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Abstract 

Wc present a new methodology to derive the positions of the Sgr B molecular clouds (MCs) along the line 
of sight, as an application study of the Galactic center diffuse X-rays (GCDX). The GCDX is composed 
of hot plasma emission of about 7 keV and 1 keV temperatures, and non-thermal continuum emission 
including the 6.4 keV line from neutral irons. The former, the Galactic center plasma emission (GCPE), is 
uniformly distributed over 1 degree in longitude, while the latter is clumpy emission produced by Thomson 
scattering and fluorescence from MCs irradiated by external X-rays (the X-ray reflection nebula emission: 
XRNE). We examined the Suzaku X-ray spectra of the GCPE and XRNE near to the Sgr B MC complex, 
and found that the spectra suffer from two different absorptions of Nn (Absl) > 10 23 H cm -2 and Ah 
(Abs2) ~ 6x 10 22 H cm -2 . Absl is proportional to the 6.4 keV-line flux, and hence is due to the MCs, while 
Abs2 is typical of interstellar absorption toward the Galactic center. Assuming that the GCPE plasma is 
spherically-extended around Sgr A* with a uniform density and the same angular distribution of the two 
temperature components, we quantitatively estimated the line-of-sight positions of the MCs from the flux 
ratio the GCPE spectrum suffered by Absl and that with no Absl. The results suggest that the Sgr B 
MCs are located at the near side of Sgr A* in the GCPE. 

Key words: Galaxy: center — ISM: molecules, clouds, structure — X-ray spectra 
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1. Introduction 

The Galactic center (GC) region is complex with many 
molecular clouds (MCs), emission nebulae, supernova 
remnants (SNR), non-thermal radio emission, and so on. 
The major component is the central molecular zone (CMZ, 
|/| < 1°; Morris & Serabyn 1996) which contains mass of 
3 x 10 7 M Q (Dahmen ct al. 1998), approximately corre- 
sponds to 10% of the whole molecular gas in the Galaxy. 
Sgr Bl and B2 are the most well-known and massive MCs 
in the CMZ. The detail of the CMZ structure has been 
studied mainly based on radio observations. Using (l,V) 
diagrams of the CO and CS lines and based on the as- 
sumption of the dynamical motions of the GC region, 
the line-of-sight distribution of MCs was investigated (e.g. 
Sofue 1995; Tsuboi et al. 1999; Nakanishi & Sofue 2003). 
On the other hand, Sawada et al. (2004) examined the spa- 
tial correlation between 2.6-mm CO line and 18-cm OH 
absorption line, and predicted the line-of-sight positions 
of MCs with no assumption of dynamical motions. 

The GC region within the CMZ is also bright in the 
diffuse X-rays called the Galactic center diffuse X-ray 
(GCDX) (e.g. Koyama et al. 2007c). The GCDX has 
been found to have thermal emission from hot plasmas, 
and hence we refer to the emission as the Galactic cen- 
ter plasma emission (GCPE). The X-ray spectrum of the 
GCPE exhibits strong K-shell emission lines from highly 
ionized atoms such as FeXXV (6.7 keV) and Sxv (2.45 
keV). Koyama et al. (2007c) found that the line-intensity 



ratios of Fe XXVI Ka (6.97 keV)/FeXXV Ka (6.7 keV) and 
FeXXV K/3 (7.88 keV)/FeXXV Ka (6.7 keV) are about 
0.33 and 0.09, respectively. These values indicate that the 
plasma is in a collisional ionization equilibrium (CIE) at 
the temperature of kT ~ 6.5 keV. On the other hand, the 
K-shell emission from lighter elements, like the Sxv Ka 
line (2.45 keV), is due to a lower temperature plasma 
of - 1 keV. The GCPE extends over the GC region of 
about 1 degree in longitude with a monotonic decrease as 
the distance from Sgr A* increases (Koyama et al. 1989, 
Yamauchi et al. 1990). 

The GCDX has another component, non-thermal emis- 
sion with the Fel Ka line of 6.4 keV from neutral irons 
(hereafter 6.4 keV-line). The diffuse 6.4 keV-line in the 
Sgr B region was firstly found by the ASCA satellite 
(Koyama et al. 1996). The Suzaku satellite deeply ob- 
served the Sgr B region and found three bright 6.4 keV 
sources, Sgr B2 (M 0.66-0.02), M 0.74-0.09 (Koyama 
et al. 2007b), and Sgr Bl (M 0.51-0.10; Nobukawa et 
al. 2008), at the locations of the MCs identified by the 
CO line and the CS line (Sofue 1995; Tsuboi et al. 1999). 
The X-ray spectra of these sources have common features: 
a large equivalent width (EW > 1 keV) for the 6.4 keV- 
line and a strong absorption (A^h > 10 23 H cm -2 ) feature 
for the continuum below ~ 4 keV. In addition, a time vari- 
ability of the 6.4 keV-line has been discovered from Sgr 
B2 (M 0.66-0.02; Koyama et al. 2008). These facts sug- 
gest that the origin of the 6.4 keV-line and the related 
continuum emission is due to fluorescence and Thomson 
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Fig. 1. XIS band images around the Sgr B region smoothed with a Gaussian kernel of a = 1'. (a) The 6.4 keV band (6.3-6.5 keV), 
which contains the Fel Ko line, tracing the dense molecular clouds, (b) Medium band (2.5-3.5 keV) showing the distribution of the 
GCDX in low energy. The CS-line contours (Tsuboi ct al. 1999) in velocity range of —50-140 km s _1 are added with green solid 
lines in (a) and (b). A hint of anti-correlation between (a) and (b) near I ~0.°75 would be due to strong absorption by the molecular 
clouds. Spectra regions are shown in black solid lines in (a); the spectrum region for M0.74-sub is the largest circle excluding the 
small circle of M 0.74-0.09. The bright point source (SAX J1748.2— 2808; Nobukawa et al. 2009) is excluded for Region 1. The scan 
profile in figure 2 is made from the large rectangle (0.4 degree X 0.1 degree) in (b). 

scattering taking place in molecular clouds irradiated by sight: a correlation study of the flux of the 6.4 keV-linc 
an external source. We refer to this component as the X- from the MCs and the line-of-sight absorption (Nb) of the 
ray reflection nebula emission (XRNE). The most prob- GCPE. 
able external X-ray source is a past outburst of Sgr A* 
(Koyama et al. 1996; Inui et al. 2009). 

In this paper, we propose a new methodology to derive 
positions of the Sgr B molecular clouds along the line of 
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Table 1. Suzaku Observations near the Sgr B Region 
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Target name 


Obs.ID 


Pointing 


direction* 


Obs.date 


Effective exposure 






a (J2000.0) 


5 (J2000.0) 


(yyyy-mm-dd) 


(ks) 


GC SgrB2 


100037060 


17 h 47 m 30?60 


-28°26'36."6 


2005-10-10 


76.6 


GC Number2 


500005010 


17 h 47 m 04?63 


-28°37'46."2 


2006-03-27 


88.4 


GC SgrB East 


501039010 


17 h 48 m 04?87 


-28°21'06."5 


2007-03-03 


96.4 


SgrD SNR 


502020010 


17 h 48 m 46?13 


-28°07'38."6 


2007-09-06 


139.1 



"The center of the XIS field of view. 



2. Observations and Data Reduction 

Deep multiple pointing observations toward the Sgr 
B region were performed with the X-ray Imaging 
Spectrometer (XIS) at the focal planes of the X-ray 
Telescope (XRT) onboard the Suzaku satellite from 
October 2005 to September 2007. The XIS system con- 
sists of three sets of front-illuminated (FI) CCD cameras 
(XIS 0, 2, and 3) and one set of back- illuminated (BI) 
CCD camera (XIS 1); each CCD chip contains 1024x1024 
pixels (1 pixel = 24 ^mx24 fim) for the 18'xl8' field of 
view. Two calibration sources of 55 Fe are installed to il- 
luminate two corners of each CCD for absolute gain tun- 
ing. The details of Suzaku, the XIS and the XRT are 
given in Mitsuda et al. (2007), Koyama et al. (2007a), and 
Scrlemitsos et al. (2007), respectively. The observations 
were made in the normal clocking mode with read-out cy- 
cle of 8 s. We screened the data by excluding the events 
observed below the low and day-Earth elevation angles at 
5° and 20°, in addition to the removal of events during 
passage of the South Atlantic Anomaly. The total effec- 
tive exposure time was about 400 ks. The observation log 
is shown in table 1. 

We performed data reduction and analysis using 
HE ADAS software version 6.5.1 and XSPEC version 
11.3.2. The calibration database was the version re- 
leased 1 on 2008-06-02. We constructed the non-X-ray 
background (NXB) data from night-Earth observations 2 , 
and subtracted the NXB from the images and the spectra. 
Since the relative gains and response functions of the FI 
CCDs (XIS 0, 2 and 3) are essentially the same, we made 
merged spectra of the FI CCDs 3 . 

3. Analysis and Results 

3.1. Overall Features of the X-Ray Images 

Figure 1 shows XIS images of the Sgr B region in the (a) 
6.3-6.5 keV and (b) 2.5-3.5 keV bands. The images were 
made by co-adding the FI and BI CCDs data, followed 
by correcting the exposure time and the vignetting effects 
of the XRTs. The data of the CCD corners illuminated 
(contaminated) by the calibration sources were excluded. 



http: / /www. astro. isas.ac.jp/suzaku/caldb/ 
http: / /www. astro. isas.ac.jp/suzaku/analysis/xis/nte/ 
The XIS 2 suddenly became unusable on 2006 November 9 pos- 
sibly due to a micro-meteoroid impact on the CCD. Therefore, 
after this epoch, the merged FI-CCD spectra were those of XIS 
and 3. 
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Fig. 2. One-dimensional profiles of the Sgr B region along 
the Galactic longitude from I = l.°0-0.°6 at b= -0.°15 ±0.°05 
(see the solid rectangle in figure lb). Gray open circles and 
black filled squares are data of the 2.5—3.5 keV band and the 
6.4 keV band (6.3—6.5 keV) with 1-cr error bars, respectively. 
The radio CS line profile (—50-140 km s" 1 ) is plotted with 
the black solid line. The flux of each data is normalized to its 
maximum value. 

To increase the visibility, we binned the images with 4x4 
pixels and smoothed them with a Gaussian kernel of a = 
1'. 

The bright clumps in figure la are SgrB2 
(M 0.66-0.02), M0. 74-0.09, and SgrBl (M 0.51-0.10). 
These are reported to be candidates of the X-ray reflec- 
tion Nebula (XRN; Koyama et al. 2007b; Nobukawa et 
al. 2008). In addition, we found a new diffuse structure in 
the 6.4 keV-line at the south of M 0.74-0.09 (hereafter, 
M0.74-sub; figure la). 

The GCDX flux in the 2.5-3.5 keV band (figure lb) is 
sensitive to absorption in Nn ranges of 10 22 -10 24 H cm -2 . 
This range corresponds to the interstellar absorption to- 
ward the GC region and the intra-cloud absorption of the 
MCs. Comparing figure la and lb, we found a hint of 
anti-correlation between the 6.3-6.5 keV (6.4 keV-line) 
and the 2.5-3.5 keV band fluxes near / = 0.°75: the 2.5- 
3.5 keV band flux is relatively weak in contrast with the 
strong emission of the 6.4 keV-line from the Sgr B MCs. 

To confirm this anti-correlation, we made one- 
dimensional flux profiles in the 2.5-3.5 keV band and the 
6.4 keV-line band along the Galactic longitude of I =1.°0- 
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0.°6 (see figure 2). For comparison, we also made profile 
and contours (in figure la and lb) of the radio CS line 
(J = 1-0) by extracting the velocity-integrated flux in the 
range of -50-140 km s" 1 (Tsuboi ct al. 1999), within 
which most of the Sgr B MCs is included. The CS-line 
flux profile is added in figure 2. For all of the profiles, 
the vertical axis represents the sum-value at the Galactic 
latitude b = — 0.°15 ± 0.°05 (see the solid rectangle in figure 
lb); the flux of each profile is normalized at the maximum 
value, and hence the relative flux range is from to 1. 

As shown in figure 2, the flux profile of the 6.4 keV-line 
shows a bump at I = 0.°65-0.°85, which is in general agree- 
ment with that of the CS line. The 2.5-3.5 keV band flux, 
on the other hand, shows an intensity drop at the bump of 
the 6.4 keV-line and the CS line, which implies absorption 
due to the MCs. Since the 2.5-3.5 keV band flux is sum of 
the GCPE behind the MCs (should be anti-correlated to 
the 6.4 keV-line flux) and the GCPE in front of the MCs 
(no correlation to the 6.4 keV-line flux), the intensity drop 
would be smeared out to only ~ 15%. Detailed study of 
the spectra is given in section 3.2. We note that the angu- 
lar size of the 6.4 keV-line bump is slightly narrower than 
that of the CS line and that of the intensity drop in the 
2.5-3.5 keV band at / - 0.°85 (section 4.3). 

3.2. Overall Features of the Full-band Spectra 

3.2.1. Complex Components in the Spectra 

We obtained X-ray spectra of the Sgr B region in the 
0.5-10.0 keV band from the 6 solid line regions given in 
figure la. These are SgrB2, M 0.74-0.09, SgrBl, M0.74- 
sub, Region 1, and Region 2. In figure 3, for brevity, we 
show 2 typical examples of the spectra. 

The spectra have two notable common features (see fig- 
ure 3). One is that the spectra exhibit many emission lines 
such as FeXXV Ka (6.7 keV), FeXXVI Ka (6.97 keV), 
Si XIII Ka (1.86 keV), Sxv Ka (2.45 keV), and Fel Ka 
(6.4 keV). These are K-shell lines from various atoms in 
various ionization states and hence the spectra are likely 
to consist of multi-temperature components, but may be 
reduced to at least three temperatures, which are respon- 
sible for the emission lines of the 6.4 keV (Fel, cool gas), 
the 2.45 keV (Sxv, hot plasma), and the 6.7 keV (FeXXV, 
very hot plasma). The other feature is that the contin- 
uum shape shows broad local minima near 1-3 keV. Since 
the absorption of the GCDX spectrum is in the range of 
6 x 10 22 -10 24 H cm~ 2 , the flux should rapidly drop below 
~ 2 keV. The presence of the local minima, therefore in- 
dicates that the full-band spectra have an additional soft 
X-ray component with no significant (small) absorption 
(the foreground emission). As a consequence, we decom- 
pose the spectra of the Sgr B regions to the following five 
components: 

(1) Very hot plasma responsible for the 6.7 keV-line. 

(2) Hot plasma responsible for the 2.45 keV-line. 

(3) Cool gas component that emits the 6.4 keV-line. 

(4) Low energy (0.5-2 keV) X-rays with small absorption. 

(5) The cosmic X-ray background (CXB). 

Koyama et al. (2007c) studied component (1) using the 



iron/nickel line flux ratios of K a and K f3 in the GC re- 
gion. They concluded that component (1) is a plasma of 
temperature kT ~ 6.5 keV in collisional ionization equilib- 
rium. The iron abundance was determined to be ~ 1 so- 
lar 4 . We therefore use an APEC model 5 for component (1) 
and express it as APEC1 with fixed temperature and 
abundances of 6.5 keV and 1 solar, respectively. 

With the Suzaku GC survey observations, many SNR 
candidates were found in the 2.45 keV-line band images. 
The plasma temperatures are nearly the same at about 1 
keV. (Nobukawa et al. 2008; Mori et al. 2008; Tsuru et al. 
2009). The abundances are slightly different in each ele- 
ment, but are consistent with 1 solar. The abundances of 
S, the most prominent Ka line emitter, are 0.6-1.3 solar 
(hereafter, the numbers in the parenthesis are in 90% er- 
ror range) for G 0.42-0.04 (Nobukawa et al. 2008), 1.2-1.5 
solar for G 359.79-0.26 (Mori et al. 2008), and 1.1-2.7 so- 
lar for G 359.41-0.12 (Tsuru et al. 2009). Removing these 
SNR candidates, we still found the 2.45 keV line prevail- 
ing in all of the GC region (component (2)). Thus we use 
an APEC model for component (2) (APEC2), assuming 
that the temperature and abundances to be ~ 1 keV (free) 
and 1 solar (fixed), respectively, the same values to those 
of SNR candidates near the GC region. 

Component (3) is emission possibly produced by the flu- 
orescence and reflection from the MCs (XRNE): the flu- 
orescent Fel Ka (6.4 keV-line) with Thomson scattering 
continuum. The spectra of XRNs in the SgrB and SgrC 
complexes were extensively studied with Suzaku (Koyama 
et al. 2007b; Nobukawa et al. 2008; Nakajima et al. 2009). 
Following these results, we use a phenomenological model: 
a power-law plus 2 Gaussian lines for the Fel Ka and Kf3 
as, 

XRNE = Ax (E/keVy F + Gaussianl + Gaussian2 

[photons cm~ 2 s _1 arcmin~ 2 ] (1) 

where the Gaussian center energy and relative intensity 
ratio of the Fel Ka, Kf3 lines are fixed to 6.4 keV and 
7.05 keV, and 1 : 0.125, according to the theoretical values 
(Kaastra & Mewe 1993). 

Component (4) is low energy (< 2 keV) emission extend- 
ing to nearly 0.5 keV, indicating small or no absorption, 
and hence the foreground emission, which is unrelated to 
the GC region. The origin is unknown but would be either 
a local Galactic plasma (Ebisawa ct al. 2008) or unresolved 
faint dM stars (Masui et al. 2009). Apart from the real 
origin, the spectra of component (4) contains faint line-like 
structure (see figure 3), we therefore used an APEC model 
(APEC3). Since the contribution of component (4) above 
~ 2 keV is small, the fitting results of the spectra above 
~ 2 keV depend very weakly on the model assumption of 
component (4). 

Component (5) is the cosmic X-ray background (CXB). 
We applied the CXB model obtained from the Suzaku 
data at the north ecliptic pole, 

4 The solar abundance in this paper is referred to Anders & 
Grevesse (1989). 

5 Astrophysical Plasma Emission Code: A model of an emission 
spectrum from collisionally-ionized diffuse gas. (Smith et al. 
2001). 
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Fig. 3. FI spectra of MO.74-0.09 (left column) and SgrB2 (right column) as examples for different models. The best-fit spectra and 
the residuals (x-distributions) of models A, B, and C are shown in the top, middle, and bottom rows, respectively. In each spectrum, 
components of the NGCE and the XRNE are shown by green and black solid lines; components of GCPE suffered by Absl X Abs2 
and Abs2 are shown by the red and blue solid lines, respectively (see section 3.2.2). 



CXB = 7.4 x 1(T 7 x (E/keVy 1Am 
[photons cm" 



(2) 



As a matter of fact, the surface brightness of the CXB 
(see also Kushino et al. 2002) in the SgrB region is far 
lower than those of the other components by more than 
two orders of magnitudes, and hence the contribution of 
the CXB is almost negligible. 



3.2.2. Model Constructions and Fitting 

For the full-band spectra of the five components, we 
tried model fitting with many free parameters. However, 
due to large errors, no significant constraint on the best-fit 
parameters was obtained for any further scientific study 
and discussion. We therefore set several constraints on 
the model based on reasonable assumptions. 

We assume that the spectral shape of the GCPE, which 
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comprises the two plasmas (APEC1 and APEC2), is the 
same in all of the Sgr B regions. We thus combine com- 
ponents (1) and (2), 



GCPE = APEC1 + APEC2 

[photons cnf 2 f 



(3) 



The temperature and abundances of APEC1 are fixed to 
6.5 keV and 1 solar, respectively. We define a parameter a, 
the normalization ratio of APEC1 and APEC2, assuming 
the common value in all regions near Sgr B (see figure 
6). Thus, free parameters for the GCPE spectra are the 
temperature and the normalization factor of APEC2, and 
a with the same value among the relevant regions. 

The photon indices (r) of the continuum of XRNs in 
Sgr B are scattered. This may be due to the limited en- 
ergy band fittings in former studies (section 4.2). On the 
other hand, the equivalent width values of the 6.4 keV-line 
are nearly the same (Koyama et al. 2007b, Nobukawa et 
al. 2008, Nakajima et al. 2009). Since the emission mech- 
anisms would be the same for all Sgr B XRNs, we assume 
that the photon index (r) and the equivalent width (EW) 
of the 6.4 keV-line are the same in all regions of Sgr B. 

Although each of components (l)-(5) should have a dif- 
ferent absorption, we truncated them into three absorp- 
tions: Absl, Abs2, and Abs3. Absi (i = 1, 2, 3) is the in- 
tegrated value along the line of X-ray path into our sight, 
which is given as exp(— Ah x <J(e))] Ah an d <J(e) ar e re- 
spectively the hydrogen column density and absorption 
cross section with solar abundance. A schematic view of 
these three absorptions is given in figure 4. 

Absl is the largest absorption (Ah > 10 23 H cm -2 ), 
which takes place in/through the MCs. Abs2 is the 
interstellar absorption to/beyond the GC of Ah ~ 6 x 
10 22 H cm" 2 (Sakano et al. 2002). Abs3 is the small- 
est absorption applied to the foreground emission (Ah < 
10 22 H cm- 2 ). 

In all 6 regions, we needed additional common compo- 
nents (4) and (5) , which are not due to the emission from 
the GC region. We combined these components and define 
them as the non GC emission (NGCE): 

NGCE = Absl x Abs2 x Abs2 x CXB + Abs3 x APEC3 

[photons cm _2 s _1 arcmin _2 ] (4) 

The term Abs2x Abs2xCXB indicates that the CXB suf- 
fers from interstellar absorption (Abs2) two times: coming 
from the front and back sides of the GCPE. Here we as- 
sume that the absorption due to the back side on the GC 
is same as that on the front side of the GC. 

Referring the schematic view of the geometry of the 
MCs (XRNE) and the GCPE, we tried two extreme cases 
(sec figure 4). (A): a MC is located at the near-side edge of 
the GCPE, then the XRNE and GCPE would have nearly 
the same absorption of AbslxAbs2. (B): a MC is located 
at the far-side edge of the GCPE; thus, the XRNE and 
the GCPE would have independent absorptions, given as 
AbslxAbs2 and Abs2. Model A and Model B, therefore, 
are given as, 

(A) : Absl x Abs2 x (XRNE + GCPE) + NGCE 

(B) : Absl x Abs2 x XRNE + Abs2 x GCPE + NGCE 



A, 



Abs2 
(Nn ~ 6) 




Foreground Emission 
Abs3 
(Ah <1) 



■ Abs2 
'(A^h ~ 6) 



Observer 



Fig. 4. Schematic view of the MCs distribution along the 
line of sight for Model C. The unit of absorption (Nn) in the 
parentheses is 10 22 H cm~ 2 . 



[photons cm 2 s 1 



arcmin 2 ] (5) 



In the fittings, the free parameters arc Absl, Abs2, 
and the normalization factors of the 6.4keV-line flux, the 
GCPE, and the APEC3. While the mixing ratio a of 
APEC1 and APEC2 (in the GCPE), Abs3, the tempera- 
ture and the abundance of APEC3, the equivalent width 
(EW) of the 6.4 keV-line and the photon index (P) for 
the power-law spectrum of XRNE are also free, but are 
assumed to be the same (common) in all relevant regions. 
We list the fitting parameters according to their properties 
in table 2. 

We performed a simultaneous fitting for the 2x6 spec- 
tra (the FI and BI spectra of 6 regions) with Model A and 
Model B. The best-fit parameters and x 2 /d.o.f are listed 
in table 3. As a results, both models were rejected by the 
X 2 /d.o.f test. Model A gave a large x 2 /d.o.f of 1.87, and 
shows large residuals in the medium range of 1-3 keV and 
the hard range of 5-7 keV (see figure 3). The best-fit val- 
ues of Absl x Abs2 (with Model A) is ~ 10 x 10 22 H cm~ 2 
(see table 3), which are close to the typical value toward 
the GC region. These values may be reasonable for the 
GCPE but could not be applied to the XRNE, because the 
absorption of each XRNE was found to be > 10 23 H cm -2 
by independent analysis of XRN (Koyama et al. 2007b; 
Nobukawa et al. 2008; Nakajima et al. 2009). 

In Model B, the apparent residuals in the medium and 
hard range found in Model A disappeared partially, and 
hence x 2 /d.o.f was improved to 1.77. We still see unnat- 
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Table 2. Summary of free and fixed parameters in the fittings for Model A, B, and C. 


Free parameters 


Absorption: 
XRNE*: 
GCPE*: 
NGCE*: 


N H (Absl), N H (Abs2) 
normalization (Fel Ka 6.4 keV line) 
normalization (APEC2), position BJ 
normalization (APEC3) 


Free parameters but common in all regions 


Absorption: 
XRNE*: 
GCPE*: 
NGCE*: 


JV H (Abs3) 

photon index T, equivalent width EWeAkcV 

kT (APEC2), a (normalization ratio of APEC1/APEC2) 

kT (APEC3), abundances Z (APEC3) 


Fixed parameters 


XRNE*: 
GCPE*: 
NGCE*: 


Fel line energy Ka = 6.4 keV, K/3 = 7.05 keV, flux ratio (K/3/Ka) = 0.125 
kT (APEC1) = 6.5 keV, Z (APEC1 & APEC2) = 1 solar 
CXB = equation 2 



* See equation 1, equation 3, equation 4 tor the definitions of XK1NF, GCPE, and NGCE, respectively, 
t For Model C only, see equation 6 and the text. 



Table 3. Results of best-fit parameters* with Model A and B. 





MODEL A 




MODEL B 




Region 


M (AbslxAbs2) 


6.4 keV-hW 


ATI (Absl) 


M (Abs2) 


6.4 keV-hW 


Sgr B2 
MO.74-0.09 
Sgr Bl 
M0.74-sub 
Region 1 
Region2 


io.7i°i 

8.4t° j 
7 4+0-1 

6 1+ - 1 

Da - 0.1 

7 7+0.1 
' ■ ' -0.3 


3.0± u; i 

1 ro+0.04 

i oc+0.02 
1 ' zd -0.02 

24+ - 01 

u - z4 -0.01 
u -^ J -0.03 


54.8+14 
40.4t 2 ° 

26.3+fg 
21.9+S'g 
18.6iJJ 
31.8±?;| 


5 8+ U ' 4 

°-°-0.1 

5 2+ 1 - 

°- z -0.8 

6.3t°;? 

s.o±8:i 


6-6^ 

4-3±8:5 

o o+O.l 
z - z -0.1 

1 63+ 03 

i -°°-0.02 

S3+ 01 

U - JJ -0.01 

7fi+ 02 

U - '"-0.04 



Region-common parameters 



kT (keV) (APEC2) 

( APECl norm \ 
" \APEC2norm'' 

Photon index T 
EW 6AkcV (keV) 
Abs3 (10 22 Hcm- 2 ) 
kT (keV) (APEC3) 
Z (solar) (APEC3) 


0.70 (0.68-0.72) 
0.47 (0.46-0.48) 
0.52 (0.47-0.59) 
1.78 (1.75-1.80) 
0.21 (0.21-0.22) 
1.02 (0.99-1.05) 
0.004 (0.002-0.006) 


0.69 (0.68-0.74) 
0.59 (0.57-0.60) 
2.16 (2.10-2.23) 
1.31 (1.28-1.34) 
0.20 (0.19-0.21) 
0.93 (0.92-0.94) 
0.008 (0.007-0.011) 




4749/2537= 1.872 


4479/2537= 1.765 



*Thc uncertainties arc at 90% confidence level. 

+ The value of AT H in the unit of 10 22 Hem" 2 . 

fFlux in the unit of 10 — 6 photons cm" 2 s~ 1 arcmin— 2 . 

§The results obtained by the simultaneous fitting of the 12 spectra (the FI and BI spectra of 6 regions). 



ural residuals between 2-4 keV. Although the x 2 /d.o.f of 
1.77 is large enough to reject Model B at more than the 
99.9% confidence level, this value is slightly smaller than 
that of Model A. At the end of this section, we discuss 
why we reject Model B from a different point of view. 

We then constructed Model C, adding a new free pa- 
rameter R (for each region), in which a MC/XRN is inside 
the GCPE with a relative position of R (0-1); the smaller 
R indicates the nearer (front) side. 

(C) : Absl x Abs2 x (XRNE + (1-R)x GCPE) 
+Abs2x (R x GCPE) + NGCE 

[photons cm -2 s _1 arcmin -2 ] (6) 

We performed a simultaneous fitting for Model C by the 



same procedure as for Model A and B. Although Model C 
has six new free parameters (i? for the six regions) more 
than Model A and Model B, the degree of freedom (d.o.f) 
of Model C (2531) is almost equal to that of Model A and 
and B (2537). On the other hand, the x 2 / d -°- f of Model C 
is largely improved to 1.34 from those of Model A (1.87) 
and Model B (1.77). 

One may argue that the spectra are contaminated by 
unresolved point sources in the Sgr B region. In fact, 
Koyama et al. 2009 estimated that the point source fluxes 
provide a significant fraction of the GCDX (GCPE + 
XRNE) emission, up to about 1/6 near Sgr A*, which 
may partially contribute to the 6.7 keV-line and the 6.4 
keV-line fluxes in the GCPE and XRNE, respectively. 
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According to the Chandra cataloged sources reported by 
Muno et al. (2006), the number of point sources per solid 
angle near Sgr B2 is ~ 10% compared to that near Sgr A*, 
while Nobukawa et al. (2008) shows that the GCDX flux 
near Sgr B2 is also ~ 10% of that near Sgr A*. We can in- 
fer that the contribution of point sources to GCDX would 
be roughly the same in the Sgr B region. The point source 
contribution may cause possible systematic errors in the 
present spectral analysis. The most serious errors, there- 
fore, should be the region-to-region fluctuation of the frac- 
tions of the point source fluxes, which may be rather small 
effects in a relative comparison of the best-fit parameters. 
At this moment, we have no information on the exact flux 
and spectrum of the integrated point sources in the Sgr B 
region. We therefore ignore the point source contribution. 

All together, taking account of possible systematic er- 
rors due to the point source contribution and those due to 
the simplified model, we regard Model C as a good approx- 
imation of the overall spectra in the Sgr B region, although 
the x 2 /d.o.f of 1.34 is statistically marginal to be accept- 
able. The best-fit common parameters for all regions are 
as follows (hereafter, the numbers in the parenthesis are 
in 90% error range). The temperature of APEC2 is 0.87 
(0.81-0.91) keV and the mixing ratio (a) in the GCPE 
is determined to be 0.27 (0.26-0.28). The photon index 
(r) and the equivalent width (EW) of the 6.4 keV-line for 
the XRNE are 1.72 (1.64-1.80) and 1.59 (1.54-1.63) keV, 
respectively. The absorption (Abs3), temperature (kT) 
and the abundance of the foreground emission (APEC3) 
are 0.17 (0.16-0.18) x 10 22 H cm" 2 , 0.85 (0.83-0.87) keV, 
and 0.011 (0.009-0.014) of solar, respectively. The other 
best-fit parameters of Model C for each region are listed 
in table 4. 

The 6.4 keV-line are likely to be fluorescence in the 
MC generated by an external X-ray source. Then the 
flux should be proportional to the product of the external 
X-ray flux and the number of iron atoms along the line-of- 
sight X-ray path in the MC, in the optically thin case at 
the energy near 6.4-7.1 6 keV (A H < 10 24 H cm" 2 ). The 
number of iron atoms is approximately proportional to 
the column density Ah of the MC. Therefore, in the case 
that the external X-ray fluxes are the same among the 
MCs, Absl should be approximately proportional to the 
6.4 keV-line flux. 

In figure 5, we plot the best-fit Model C values of Absl 
and Abs2, as a function of the 6.4 keV-line flux. For a com- 
parison, we also plot the results of Model B. A notable fact 
to support Model C is that Absl is in good proportional- 
ity to the 6.4 keV-flux. This feature is in contrast to the 
results of Model B, where the proportionality is not very 
good in the low-flux range. Most importantly, Model B 
requires an extra off-set of about 2 x 10 23 H cm~ 2 , which 
is un-realistic. This leads us to exclude Model B in ad- 
dition to the x 2 /d.o.f test. We estimated Ah values of 
the Sgr B MCs from the data of the CO line (Sofue 1995; 
Dahmen et al. 1998) and the CS line (Tsuboi et al. 1999). 
They arc in the range of Ah = 10 23 -10 24 H cm~ 2 , in good 
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E = 7.1 keV is the K absorption-edge of neutral iron. 



1 2 3 4 5 6 7 

6.4 keV-line flux [10~ 6 photons cm" 2 s" 1 arcmiri 2 ] 

Fig. 5. Correlations between the 6.4keV-line flux and ab- 
sorptions (made from data listed in table 3 and table 4). 
Results of Model C and Model B are plotted by the black and 
gray marks; filled squares and open circles represent Absl and 
Abs2, respectively. 

agreement with the Model C prediction. 

As shown in figure 5 for Model C, we found that Abs2 is 
almost constant at Ah ~6 x 10 22 H cm~ 2 . This is consis- 
tent with the view that Abs2 is interstellar absorption (see 
figure 4) toward the GC of about Ah - 6 x 10 22 H cm -2 
(Rieke ct al. 1989; Sakano ct al. 2002). 

4. Discussions 

4-1. Justifications of Model Constraints & Assumptions 

In the previous section, we introduced a complicated 
model to described the spectra of the GCDX. In order 
to obtain meaningful values (no large errors) of physi- 
cal parameters, we further added some constraints and 
assumptions. If, for example, we let the normalization 
ratio a (APEC1/APEC2) to be free for all regions, we 
can not give a significant constraint on the absorptions 
of GCPE and/or the temperature of APEC1, because the 
two plasma emission become nearly the same in the band 
near 2 keV. We summarize the adopted constraints and 
assumptions below for further discussions: 

• The 6.5 keV plasma and ~ 1 keV plasma, which 
comprise the GCPE, have nearly the same angular 
distribution near the Sgr B region, (see figure 6) 

• The abundances of the GCPE are nearly 1 solar, the 
same value as those of nearby new SNR candidates. 

• The XRNE has a common photon index (r) and 
EWe. 4kcV near the Sgr B region. 

• The point sources contribution to the GCDX is small 
and hence could be ignored in the Sgr B region. 

The most fundamental constraint is the first item; the 
factor a is constant from region to region. This implies 
that the angular distributions of the very hot (6.5 keV) 
and hot (~1 keV) plasmas are nearly the same near Sgr B. 
In order to judge this constraint, we show the flux distri- 
bution of the 6.7 keV and 2.45 keV lines along the galactic 
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Region 


(Absl) 


(Abs2) 


6.4 kcV-linc* 


APEC3 norm§ 


GCPE norm" 


Position R# 


Sgr B2 
MO.74-0.09 
Sgr Bl 
M0.74-sub 
Region 1 
Rcgion2 


36.5™ 
22.9+^ 
12-3±iig 

7 O+0.2 
' - z -0.2 

19 n +0 - 6 

lz - u -0.8 


7 o +0 - 2 

' - u -0.3 
c o+O.l 
O.O_ 2 
o fi +0.2 

6 3+ 01 
°-°-o.i 


r -,+0.3 
°- x -0.3 
q q+0.2 
O.O_ 2 

i.8±g:i 

u - ZD -o.oi 
n 47+0.03 

u -*' -0.03 


09 +u ui 

u - uy -0.01 

n n+0- 01 
o.n±S:8i 

O.OQOtfe 
u.iuo_ 004 

u - iou -0.036 


4 n+0.3b 

2 67+0- 29 

z -°' -0.28 

3 4fi+°- 25 
°-^°-0.38 

1 66+ 010 
^""-o.oq 

1 4n+°- 05 

i - 4u -0.04 

9 ss+o- 11 

Z - SS -0.18 


44+" 0b 

u -^-0.06 

n 3i + 04 

n 41+0-06 

u -^ i -0.06 

26+ 04 

u - ZD -0.04 

n 1 5+ 001 

34+0- 04 

u -°*-0.04 


xVd.o.F* 






3382/2531= 1.336 







t The value of jV H in the unit of 10 22 Hem" 2 . 

i Flux in the unit of 10 -6 photons cm - 2 s — 1 arcmin -2 . 

§ The normalization factor for the foreground emission (c.f APEC3 in equation 4) of the APEC model re- 
normalized with the region size A [arcmin 2 ]; it is expressed as 7.07 X 10 — 11 /(AttD 2 A) EM [cm - 5 arcmin — 2 ], 
where D and EM are the distance to the source [cm], and the emission measure [cm -3 ], respectively. 
II In same expression as §, but the normalization factor of APEC2 in the GCPE (c.f equation 3 and text). 
* Ft, is the fraction of the GCPE suffered by Abs2, which indicates the positions of MCs (see equation 6 and 
figure 7). 

** The results obtained by the simultaneous fitting of the 12 spectra (the FI and BI spectra of 6 regions). 



longitude of I = l.°2-0.°4 at b = -0P15 ±0.°05 (figure 6). 
The 6.7 keV and 2.45 keV lines are strong characteristic 
lines that well represent, respectively, the distributions of 
the 6.5 keV plasma and the ~ 1 keV plasma with the least 
systematic errors. 



1 



6.7 keV Line 
2.45 keV Line 




0.01 



1.2 1.1 1.0 0.9 0.8 0.7 0.6 
Galactic longitude [degree] 



0.5 0.4 



Fig. 6. Flux distribution of the 6.7 keV line (6.5 keV plasma; 
dashed line) and the 2.45 keV line (~ 1 keV plasma; solid 
line) along the galactic longitude from I = l.°2-0.°4 at 
b= — 0.°15 ±0.°05. Two parallel dotted lines guide eyes to 
see the same distribution of the 6.5 keV plasma and the ~ 1 
keV plasma. 

As is found in figure 6, the 2.45 keV line and the 6.7 keV 
line have almost the same slope (angular distribution). 
The small drop of the 2.45 keV-line flux near at I = 0.°85- 
0.°6 would be due to absorption by the MO.74-0.09 MC 
and the drop at I = 0.°6-0.°4 is due to the absorption by 
the Sgr Bl MC. This figure, therefore, indicates that the 
emission of the ~1 keV plasma (APEC2) and the 6.5 keV 
plasma (APEC1) have a common normalization ratio (a) 



for the GCPE near the Sgr B regions. The justification 
of the other constraints and assumptions listed in the last 
three items are separately discussed in the previous sec- 
tions. We thus conclude that Model C with the above 
listed constrains and assumptions provides a practical ap- 
proximation for the overall spectra in the Sgr B region, 
although the x 2 /d.o.f of 1.34 is statistically marginal to 
be acceptable. 

4-2. Line-of-Sight Positions of MCs 

From the results of section 3, we conclude that the MCs 
of Sgr B should be located in the GCPE (Model C), nei- 
ther in the near-side edge (Model A) nor in the far-side 
edge (Model B) of the GCPE. For simplicity, we assume 
that the volume emissivity of the GCPE is approximately 
a uniform sphere with a uniform plasma density on the 
Galactic plane. Koyama et al. (1989) and Yamauchi et 
al. (1990) reported that the 6.7 keV line distributed with 
a Gaussian shape of 1.8 deg (FWHM) (in parallel to the 
Galactic plane). We therefore take the boundary (radius) 
at \l\ =1.25 deg, within which 90% of the GCPE is in- 
cluded. The center of the GCPE is approximately Sgr 
A* (at I = -0.°056; Yusef-Zadeh et al. 1999). Then the 
parameter R gives the line-of-sight positions of the MCs 
within the GCPE 7 , where R = and 0.5 indicate the ante- 
rior border line and the center line (Sgr A* is on this line) 
of the GCPE, repectively. We plot the best-fit position 
parameter (R) for the 6 regions (indicated in figure la) in 
the schematic view of figure 7. 

As demonstrated in figure 7, the MCs in the Sgr B re- 
gion are generally on the near-side with respect to Sgr A*. 



The emission measure EM can be written as J n c rindV = 

UeUfiS J dl oc n e nnSL, where L is the length of the GCPE 
in the line of the sight, S and n e nn are the projected area and 
the product of densities of electrons and hydrogens, respectively. 
Then the parameter R is the fraction of £Afg°p E /EM^ pE , 
which is equal to ^gcpe/^GCPE : ^ e rela tive line-of-sight po- 
sition of MCs in the GCPE. 
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[degree] 
0.75 0.5 




*1 

Observer 



Fig. 7. Relative positions of each MC (filled circles) along 
the line of sight, where the horizontal lines of R = 0.0 and 0.5 
imply the anterior border line and the GC line (Sgr A* is on 
this line), respectively. The circles filled with darker colors 
indicate MCs of the larger JVh values. 

Also if all of the MCs/6.4 kcV clumps, Sgr Bl, Sgr B2, 
MO. 74-0. 09, and M0.74-sub are the constituent parts of 
the CMZ, the results indicate that the CMZ at the pos- 
itive Galactic longitude has a bar-like distribution that 
generally inclines toward us. Although no further quanti- 
tative comparison is available with the present data, the 
inclining trend of the Sgr B MCs toward us is consistent 
with the radio molecular observations results based on the 
(l,V) diagrams study by Sofue (1995) and the spatial cor- 
relation study between 2.6-mm CO line and 18-cm OH 
absorption line by Sawada et al. (2004). 

It may be suspicious that the best-fit values of R in 
Model C indicate that the MCs are all closer to the near- 
side edge of the GCPE rather than to the far-side edge, 
while x 2 /d.o.f. = 1.872 of Model A (MCs at the near-side 
edge) is rather worse than x 2 /d.o.f. = 1.765 of Model B 
(MCs at far-side edge). As shown in equation 5, Model 
B have the extra-part Abs2 x GCPE which Model A does 
not have. This gives a slightly better x 2 /d.o.f. of Model B 
(1.765) than Model A (1.872). On the other hand, the 
improvement of Model C (x 2 /d.o.f.= 1.336) from Models 
A (and B) is far larger than that between Model A and 
B. 

4.3. A Unified Picture of the XRN 

According to former studies (Koyama et al. 2007b; 
Nobukawa et al. 2008), the equivalent width of the 6.4 
keV-line and the continuum photon index (EW, V) of the 
Sgr B XRNs have been reported to be (1.1 keV, 3.2) for 



Sgr B2, (1.55 keV, 1.4) for MO.74-0.09, and (1.4 keV, 1.8) 
for Sgr Bl. The T values are significantly scattered among 
the XRNs, but the variation of EW is not very large. 
The GCPE flux is position-dependent, which means that 
a proper background-selection is difficult. The large scat- 
ter of r values in the former studies may be due to the 
nearby GCPE-subtraction and would also be due to the 
band limited (e.g. 5-8 keV) spectral fitting. In this work 
(section 3.2 with Model C), the EW and T values are esti- 
mated by an independent method, a simultaneous fitting 
of a wide energy band (0.5-10.0 keV), in which nearby 
GCPE was not subtracted but directly included as the 
GCPE model. We also note that these uncertainties do 
not largely affect the EW values. In fact, We found that 
the averaged EW ~ 1.6 keV, which is consistent with the 
previous studies of the Sgr B XRNs. Koyama et al. (2009) 
studied the GCDX near the Sgr A region < 0.°2), and 
found that the T of the continuum flux is 1.9; this re- 
duced to 1.4 if we exclude the point sources contribution. 
Nakajima et al. (2009) studied the 6.4 keV clumps in the 
Sgr C region and found T to be in range of 1.6-1.9. We 
determined the common T in Sgr B to be ~ 1.7; this value 
agrees with the results in Sgr A, B, and C. Thus, we in- 
fer that the XRNE may have a unified photon index of 
1.4-1.9 in the overall GC region. 
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6.4 keV-line flux [10~ 6 photons cm" 2 s" arcmin ] 



Fig. 8. The best-fit absorption (Absl) vs the distance-cor- 
rected 6.4 keV-line flux from results of Model C in table 4 
(open triangles), where the correction factors are estimated 
from the distances between MCs and Sgr A* with figure 7, 
and are normalized to that of Sgr Bl. 

As discussed by Koyama et al. (2008) and Nobukawa 
et al. (2008), a possible external X-ray source is Sgr A*, 
therefore the X-ray flux at the Sgr B MCs should be 
proportional to LqD~ 2 , where Lq and D are, respectively, 
the luminosity of Sgr A* and the distance between 
Sgr A* and MCs. Thus the flux of the fluorescent 6.4 
keV-line should be modified by this distance effect. In 
figure 8, we plotted distance-corrected 6.4-keV fluxes as 
a function of Absl, where the distances are estimated 
using figure 7. We then fitted the correlation between 
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the best-fit absorption (Absl) and the 6.4 keV-line flux 
(distance-corrected), and obtained, 

A H (Absl) = 3.99 (±0.56) x / 6 . 4k oV [10 22 H cm" 2 ], 

^6.4kcV : [10 -6 photons s" 1 cm" 2 arcmin -2 ] (7) 

This relation is generally in agreement with the fluorescent 
process due to irradiation by Sgr A* with the past lumi- 
nosity of 2-3 xlO 39 ergs s" 1 . If La was time constant, the 
relation between the distance-corrected 6.4 keV-line flux 
vs. Absl should be in good proportionality. However, we 
see some data scatters from the proportionality relation 
(figure 8). It can be explained by assuming a time variable 
flux of the irradiating source Sgr A* (Koyama et al. 2008; 
Inui et al. 2009). 

4-4- The Correlation of CS line and 6.4 keV-line Flux 

As shown in figure 2, the 2.5-3.5 keV band flux shows 
an intensity drop at the bump of the 6.4 keV-line and 
the CS line. The anti-correlation to the 2.5-3.5 keV 
band flux is clearer in the CS line than that in the 6.4 
keV-line. This may be because that the CS line traces all 
(GC + foreground) MCs in the line of sight from which 
the GCPE may suffer absorption, while 6.4 keV-line only 
traces MCs near the irradiating source in the GC region. 
The better anti-correlation of the CS line against the 
2.5-3.5 keV band flux indicates that a large amount of 
MCs, identified by the velocity-integrated CS lines but 
not bright in 6.4 keV-line, is located in the near side 
of the GCPE. The angular size of the 6.4keV bump is 
slightly narrower than that of the CS line (I ~ 0.°85 and 
0.°63). The narrower angular size of the 6.4 keV bump 
than that of the CS line may also be explained if Sgr 
A* is the external X-ray source to produce the 6.4 keV 
clumps. The surface brightness of fluorescent X-rays is 
faint at the near and far sides of the cloud-limb because 
of small Ah values, but becomes brighter when the 
external X-rays go into the large Ah (> 10 23 H cm" 2 ) 
region in the MCs. Another possible reason of the faint 
6.4 keV-line flux at the far-side cloud limb from Sgr A* 
is due to the time variability of the Sgr A* X-rays; an 
X-ray front of possible X-ray outbursts about 300 years 
ago (Koyama ct al. 2008; Inui et al. 2009) does not yet 
reach at the far-side limb of I ~ 0.°85. 
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